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ABSTRACT 

We discuss the possibility to observe the products of dark matter annihilation that 
was going on in the early Universe. Of all the particles that could be generated by this 
process we consider only photons, as they are both uncharged and easily detectable. 
The earlier the Universe was, the higher the dark matter concentration n and the 
annihilation rate (proportional to n 2 ) were. However, the emission from the very early 
Universe cannot reach us because of the opacity. The main part of the signal was 
generated at the moment the Universe had just become transparent for the photons 
produced by the annihilation. Thus, the dark matter annihilation in the early Universe 
should have created a sort of relic emission. We obtain its flux and the spectrum. 

If weakly interacting massive particles (WIMPs) constitute dark matter, it is 
shown that we may expect an extragalactic gamma-ray signal in the energy range 
0.5 - 20 MeV with a maximum near 8 MeV. We show that an experimentally observed 
excess in the gamma-ray background at 0.5 - 20 MeV could be created by the relic 
WIMPs annihilation only if the dark matter structures in the universe had appeared 
before the universe became transparent for the annihilation products (z ~ 300). We 
discuss in more detail physical conditions whereby this interpretation could be possi- 
ble. 
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1 INTRODUCTION 

Though the cosmological measurements (Tegmark ct al. 
120041 ') show that there must be approximately five times as 
much dark matter as all baryon one, its physical nature re- 
mains unknown. Now the most commonly used hypothesis 
is that it consists of some elementary particles generated in 
the early Universe (hereafter we will call them Dark Matter 
Particles, DMPs). These particles are uncharged and do not 
interact strongly; there are telling arguments to believe that 
they were cold (v <C c) in the epoch when the relic radiation 
was generated. It is worth mentioning that a particle with 
suitable properties hasn't been discovered yet, in spite of no 
lack of theoretical candidates predicted by various quantum 
field theory models. 

If the premise is true, the dark matter is a mixture 
of equal quantities of particles and antiparticles, and they 
must collide and annihilate wherever the dark matter is 
present. Experimental observation of such a process would 
give us some valuable information about the DMP nature. 
In this article, we consider the dark matter annihilation in 
the epoch near the relic radiation formation (z ~ 1000). 



At that time, the average dark matter particle concentra- 
tion n oc (z + l) 3 was nine orders higher than now and four 
orders higher than in our Galaxy near the Sun system. So 
we may expect that the annihilation, the rate of which is 
proportional to n 2 , was very intensive in that epoch. 

Of all the particles that can be generated by the dark 
matter annihilation we will consider only photons, as they 
are both uncharged and easily detectable. Uncharged parti- 
cles do not interact with the magnetic field of the Galaxy, 
which allows one to measure the extragalactic background 
reliably enough. 

The earlier the Universe was, the higher the DM den- 
sity and the annihilation rate were. However, the emission 
from the very early Universe cannot reach us because of the 
opacity. The main part of the signal was generated at the 
moment the Universe had just become transparent for the 
photons produced by the annihilation. Later the dark matter 
density rapidly dropped, decreasing the signal. The moment 
(and its redshift) depends, of course, on the characteristic 
energy of the photons, in other words, on the nature of the 
dark matter. 

Thus, the DM annihilation in the early Universe should 
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have produced a sort of relic emission. We obtain its flux 
and the spectrum. On the one hand, such an emission can 
be detected in the spectrum of the extragalactic background. 
A distinguishing feature of this radiation should be its high 
isotropy. On the other hand, the absence of such a signal 
can impose a severe limitation on a dark matter model. 

In calculations we assume that the dark matter is ho- 
mogeneous. This is the simplest, but by no means the most 
natural supposition. Possible amplitude of the dark matter 
density perturbations and their influence on the annihilation 
signal will be discussed at the end of the article. 

In the second part of the article, the case of the most 
popular dark matter candidate - Weakly Interacting Mas- 
sive Particles (WIMPs) - is considered in more detail. We 
demonstrate that the relic signal from the WIMPs annihila- 
tion might have been already observed. 



2 CALCULATIONS 

The metric of a homogeneous isotropic universe can be rep- 
resented as ds 2 — c 2 dt 2 — a 2 (t)dl 2 , where dl is an element o f 
three-dimensional length l)Landau fc Lifshit3 lanv edition!) . 
In ag reement with the observations jGorbunov fc Rubakovl 
2008) we assume zero three-dimensional curvature of our 
Universe. We introduce polar coordinates with the centre at 
the local observer. Then the metric can be written as: 

ds 2 = cdt 2 - a 2 (t) [dr 2 + r 2 (d(, 2 + sin 2 C df) ] (1) 

We choose the normalization of a so that the contemporary 
value ao = 1 (hereafter the subscript is used to describe the 
present-day values of quantities). Then z and a are related 
by the following equation: 



z + 1 = - 

a 



(2) 



The number of annihilations in the volume dV in an interval 
of the proper time dr is (we should remind once again that 
we consider a homogeneous dark matter distribution) 



- (uv) n dr dV 



(3) 



Here the multiplier = takes into account that the annihila- 
tion is only possible if a particle collides with an antiparti- 

Here we should make two temporary simplifying suppo- 
sitions: we assume that one act of annihilation produces one 
photon of fixed energy (3 and that the Universe is transpar- 
ent for the photons (later we will generalize the result taking 
into account opacity and arbitrary spectrum of the emitting 
photons). Let us consider a three-dimensional space element 
that is part of a spherical layer [r; r + dr] viewed by the lo- 
cal observer from a solid angle do. If the effective area of the 
detector is dS, the number of photons arriving at it in an 
interval of time dt is: 



P= M 4dtdodrdS 
8tt a 2 



(4) 



We do not consider the situation when the dark matter particle 
is identical to its antiparticle. In this case the multiplier should 
be 1 instead of ^ 



Photons emitted with the fixed energy (3 arrive at the ob- 
server with a smaller energy e in consequence of the redshift; 
f3 and e are related by the usual law 

£ 



(5) 



where a corresponds to the moment of emission. The bigger 
r is, the smaller a was. In order to find the relationship 
between them, we should write the equation of motion of a 
photon, that is ds — 0. From ^ we obtain for the radial 
motion cdt — adr; whence it follows: 



c , c dt c , 

dr = — dt — — — da = — - da 
a da 



a 



an 



(6) 



The expansion of the Unive rse after the radiation-domina ted 
epoch can be described as (|Gorbunov fc Rubakovll2008l ): 



H \ n 



(7) 



As we will see, the main part of emission appears when - 3> 
10, and we may neglect £7a in this formula. Thus, 



Substituting © and © into (QJ we get: 

P = {aV) ^4^dtdodadS 

or, if we replace a with e/(3, in accordance with |(3J, 
5 

^ dt do d^- dS — 
dt do de dS 



(8) 



(9) 




(10) 



So we obtain the sought-for spectral intensity Q of the pho- 
ton flux, i.e., the number of photons that come to the local 
observer from unit solid angle per unit time per unit of area 
per unit energy interval is: 



Q 



(<JV) UqC f3yff3 



(11) 



One can see that the integral of this equation over the energy 
diverges when e — > 0. Of course, to accomplish this, arbitrary 
early Universe should be transparent for the photons, which 
is not the case. The photons interact with the baryonic mat- 
ter; let us assume that this process can be characterized by 
some averaged cross-section N. Then the number of inter- 
actions of the primary photons per an element of physical 
length dl = adr is 



— - = —'&n b dl = — H — 7t adr = 
N a 3 



-Kcn b ^ (12) 
a A a 

Here n b is the baryon concentration n b = no/a 3 ; we also used 
([6|. In order to integrate this equation, we suppose that H 
is constant. This assumption is usually not very rigorous. 
Actually, it means that we neglect the cross-section depen- 
dence on the photon energy (which strongly changes during 
its propagation via the redshift) and possible phase transi- 
tions in the Universe. However, since the baryon concentra- 
tion n b rapidly decreases with the Universe expansion, only 
the moments just after the emission contribute noticeably to 
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the opacity. So for the constant cross-section approximation 
to be feasible we need only that the dependence of the cross- 
section on the energy is not very abrupt, which is usually 
the case. The only phase transition which occurs in the Uni- 
verse near z ~ 1000 is recombination. It changes drastically 
the cross-section of photons with an energy < 10 eV, the 
approximation K = const is obviously wrong in this case. 
However, the energy of photons generated by the dark mat- 
ter annihilation is most likely considerably higher, and the 
recombination does not affect much their cross-section. 

Here we also suppose that the photon interaction 
with the baryonic matter is pure absorption and inter- 
acting photons just disappear. This assumption is more 
rough: the contribution of s cattering can be significant (see 
IZdziarski fc Svensso 2 [|l989h for details), and we underesti- 
mate the soft part of the spectrum. However, the main char- 
acteristics of the spectrum (the characteristic energy and the 
intensity) are determined by the moment the Universe be- 
comes transparent for the photons produced by the annihi- 
lation, and this moment does not depend much on whether 
it is absorption or scattering. Consequently, the assumption 
is quite acceptable for an estimation model. Below we will 
improve it. 

Moreover, we introduce a new variable \ for the scaling 
factor of the Universe. It has exactly the same meaning as 
o, but a represents the scaling factor at the moment when 
the photon is emitted, while x represents the scaling factor 
changing during the photon propagation to the observer. 
Consequently, \ changes from a to ao = 1. Then (|12[) can 
be rewritten as 



dN b d X 



-Hcn|) dx 



(13) 



X 3 X X 2 VxHoV^ n + riAX s 

We used (0 to obtain it. Since the right part of this equation 
rapidly drops with x growing, we can simplify it: 



dN b 

— — = —Hcn - „ 

N X 2 VxHoV^ 

Now equation (|12p can be easily integrated. 



exp 



d X 



N initial V " J a X^ ' \/X HoV^r, 

It is convenient to introduce a new constant 

Kcno 



From (|15[) we obtain: 

N 



Nt, 



— exp 



P 



2 P } 2 

3 a^/a 3 X\fx 



(14) 



(15) 



(16) 



(17) 



As it has been already mentioned, x changes from a to ao = 
1. Since a < ao = 1, we may omit the last term of fjlTj) . 
Replacing a with e/ (5 according to ([5]) we obtain 



N 



N initial 



exp 



3 Eyfe ) 



(18) 



In order to allow for the opacity of the early Universe, we 
should multiply equation (jTTJ by this exponential factor. 
Finally, we obtain: 



{av)n 2 c Pyffl ( 2 /V^\ 



(19) 



Since we supposed that one act of annihilation produces one 
photon of fixed energy f3, this equation actually specifies 
only a Green's function. If the produced photons have some 
distribution f(/3)d/3, we should convolute Q(s,j3), defined 
by (|19[) . with the distribution / to obtain the spectrum that 
appears to the viewer: 



/ 



0(e) = / Q{e,0)f(j3)d/3 



(20) 



Nevertheless, the spectrum given by (|19[) is fairly wide, and 
provided that the distribution /(/?) is not too broad the 
shape of the resulting spectrum remains similar to (|19p . 



3 THE CASE OF WIMP DARK MATTER 

Equations (|19p . (I20|l are applicable in quite general cases, 
but they were obtained on a rather rough assumption that 
the photon interaction with the baryonic matter is pure 
absorption. In order to make the model more realistic, we 
should consider a certain DMP candidate. If we make a very 
natural supposition that the DMPs were in thermal equilib- 
rium with other particles in the early Universe, we can esti- 
mate their annihilation cross section l|Kolb fc Turner! 1 199d : 
Ide Boer et ai1l2005T l: 



(av) 



2 • 10~ 27 (cm 3 /s) 

O.DAlh 2 



(21) 



For the present value (|Spergel et al.l l2003h Of D.DMh = 
0.113 we obtain (av) ~ 2 • 10~ 26 (cm 3 /s), i.e., a cross sec- 
tion typical for weak interactions. Besides, the DMP must 
be massive for the dark matter to be cold in the epoch when 
the relic radiation was generated. These are the telling rea- 
sons to believe that the dark matter consists of Weakly In- 
teracting Massive Particles (WIMPs). They are expected 
to annihilate into fermion-antifermion or gauge boson pairs 
with a large fraction of quark-antiquark pairs. A WIMP 
pair annihilation finally leads to 30 — 40 photons gener- 
ation (in the fragmentation process, mainly from 7r° de- 
cays). A greater p art of the photons has energy in the range 
from 2 to 4 GeV ijjungman. Kamionkowski. fc Griestlll996l ; 
Ide Boer et alj|2005h . 

Photon propagation at cosmological redsh i fts wa s ex- 
tensively investigated in Zdz iarski fc Svensso n| (|l989h . For 
the photons of energy 2-4 GeV the main channel of interac- 
tion with the baryonic matter is pair production on atoms 
and ionized matter (jA — ► Ae + e~). As we shall see be- 
low, the Universe becomes transparent to 2 - 4 GeV pho- 
tons remarkably after the recombination. The pair produc- 
tion cross-section per one baryon averaged over the chemical 
composition (mass fraction of hydrogen and helium are 75% 
and 25%, respectively) is 



N = 4.64 are In 



513^ 
/i + 825 



(22) 



Here \x 



As we can see, the cross- 



section depends only slightly on the energy of photon. 
N = 1.63 • 10~ 2e cm 2 if ju = 4000. Now we should substi- 
tute the obtained numerical values for {av) and H into (|19p . 
but before we can improve the model. To do this, we roughly 
estimate the number and the energy of photons generated by 
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the process fA — > Ae + e~ . It produces one electron and one 
positron: each of them carries away approximately half of the 
energy of the primary photon. Subsequently, the positron 
annihilates with an electron generating two photons. Hence, 
their energy is, on average, a quarter of the primary one. Of 
course, these photons should have rather a wide distribution 
but we neglect it and think of that one primary photon of en- 
ergy e produces two secondary photons of energy e/4. Since 
the cross-section only slightly depends on the photon 
energy, we suppose that a secondary photon interacts with 
the matter with the same cross-section H as the primary one. 

We consider the secondary photon scattering as pure 
absorption. Indeed, the photons generated by them via the 
pair production have energies ~ 100 MeV. In this energy 
range the Co mpton scattering becomes the main process 
of interaction (|Zdziarski fc Svensson|[l989l ) . Its cross-section 
rapidly grows with the photon energy decreasing, and low- 
energy photons are scattered time and again, quickly loos- 
ing the energy. Thus, a scattering of a photon of energy 
^ 400 MeV is actually equivalent to its absorption. 

The number of the secondary photons is symbolized by 
A^2 • In order to obtain a differential equation describing their 
propagation we can use the same procedure as in (|13H19p . 
The relationship (|14p acquires the form 



dN 2 = 2Np ■ 



N 2 p- 



(23) 



x 2 Vx x 2 Vx 

The second term here describes the secondary photon ab- 
sorption by the substance, the first one represents their pro- 
duction via the above-mentioned mechanism (the multiplier 
2 takes into account that one act of a primary photon ab- 
sorption produces two secondary photons). The number of 
primary photons N is given by (I17|l . As the boundary con- 
dition we should use the fact that there were no secondary 
photons when the primary had been just generated 



N 2 



= 



(24) 



The solution of 



N 2 — Ni nit i a l — j 



x exp 



satisfying the conditions is 
1 1 

2 P 2 p 



(25) 



^xVx 3a^/a / 

In order to calculate the value for the observer we should 
substitute x = a o = 1. Since a«l 



N 2 



4 1 / 2 p \ 

o a^j a \ 6 a\J a J 



(26) 



Here Nmiual is the initial number of primary photons. As 
in the case of the derivation of equation (I19|l . we substitute 
© for N xnitia i 



P 2 = -Mlfc Jiexp 



3 a^fa 



dtdodadS (27) 



We suppose that the energy of a secondary photon is a quar- 
ter of the energy of the primary one. To account for this 
4e 

effect, — should be substituted for a instead of Finally, 
we obtain for the flux Q 2 of secondary photons: 



0: 



{av)nlcp p 3 ( p PVP\ 



(28) 




E [GeV] 

Figure 1. The spectra of photon s generated by variou s WIMPs 
annihilation channels, taken from lde Boer et al.l (I2005T ) (the nor- 
malization is arbitrary, the WIMP mass is 100 GeV) . The bb chan - 
nel is most important ( Jungman, Kamionkowski, & Griestlll996t) . 



The total signal is the sum Q + Q 2 . If the annihilation gener- 
ates photons with some distribution / (one act of annihila- 
tion produces f(/3) dp photons in the energy interval d/3), we 
should convolute the sum with the distribution, by analogy 
with (f2H)) : 



Q{e) = 



j (Q(e,P) 



+ Qa(e, /?))•/(/?) dp 



(29) 



Figure [1] (taken from Ide Boer et ail (|2005T )') represents 
the photon spectra caused by various WIMPs annihilation 
channels. In the interval 1 — 5 GeV the upper curve (which 
describes the main bb channel) is well approximated if the 
distribution function f(P) oc /3 _1 exp(— 0.15/3). We shall use 
the following / (P is expressed in GeVs): 



/(/?) 



26.7-/?- 1 exp(-0.15/3), P 6 (1 - 5) GeV 
0, P £ (1 - 5) GeV 



(30) 



We normalized it considering that one act of annihilation 
generated on average 30 photons. 

Strictly speaking, the annihilation cross-section and the 
relative contribution of various channels depend on the en- 
ergy of the DMPs collision, which can be importa nt, if the 
annihilation occurs, for instance, near a black hole (jBaushevi 
2009). However, in the considered case the relative velocity 
of the DMPs is very small, and we may neglect this effect. 



4 COMPARISON WITH OBSERVATIONS AND 
DISCUSSION 

In order to fulfil the calculations, the cosmological pa- 
ram eters should be concre t ized. We use the following set 
fsee iGorbunov fc Rubakovl |2008) and references therein): 
Ov = 0.75, Q m = 0. D M + fib = 0.25, n b = 0.042 (of 
course, 57a + Qdm + = 1), the Hubble constant Ho = 
2.4 ■ 10 _18 s _1 , the relic radiation temperature 2.725 K, the 
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Figure 2. The spectrum of the relic gamma-ray background, cal- 
culated according to 129I I, II30I I. 1191 . and 128H . Separate contri- 
butions of the primary and the scattered photons are represented 
by the dashed and the dotted lines, respectively. 



baryon-photon ratio r\ = nt/nph — 6.1 • 10~ 10 . We obtain 
the present baryon concentration rig = 2.5 ■ 10 -7 cm -3 . The 
DMP concentration, with DMP mass taken as Mdmp = 
50 GeV, is n = 2.5 ■ 10~ 8 cm~ 3 . Equation (gTJ gives 
(<rv) ~2 - 1(T 26 (cm 3 /s). 

The influence of the relic annihilation on the ionization 
history of the universe is negligible. Indeed, for the recombi- 
nation epoch (z ~ 1200) we have: n b ~ rigz 3 = 4.3-10 2 cm -3 , 
n ~ noz 3 — 43 cm -3 . The number of annihilations in a unit 

volume per unit time is -^{cv) n 2 — 1.9 ■ 10 _23 cm _3 s _1 . 

In the characteristic time (of the order of the hydrogen ion 
recombination time at that epoch 435 years ~ 1.4 • 10 10 s 
|Gorbunov fc Rubakovll2008l )) we have 2.6 ■ 10~ 13 annihila- 
tions that produce 2.6- 10 -2 eV of energy. A hydrogen atom 
ionization requires ~ 14 eV. So, even all the energy pro- 
duced by the annihilation is enough to ionize only 4 ■ 10 -4 % 
of atoms. 

We take N = K(2GeV) = 1.63- 10~ 26 cm 2 J22J|. Then the 
constant p is equal to 1.0 • 10~ 4 . The resulting spectrum of 
the relic gamma-ray background, obtained with the aid of 
J25J|, (133), OH, and {2"g]l. is represented in Figure [2] Instead 
of the photon flux Q we have plotted the quantity e 2 * flux 
(that is e 2 ■ Q), which is traditionally used in experimen- 
tal data picturing. The primary (I19|) and the secondary Q280 
photon separate contributions are represented by the dashed 
and the dotted lines, respectively. The contribution of sec- 
ondary photons is small in these coordinates, even though 
their total number is very large. Our main conclusions are 
that the spectrum grows up to ~ 8 MeV, and the bulk of the 
signal lies in the range from 0.5 to 20 MeV. Characteristic 
redshift of the relic gamma-rays can be easily calculated. 
According to (|19l) . the quantity s 2 ■ Q has its maximum 
at 4 = a = (2p) 2 / 3 . According to it corresponds to 
z ~ 300. 

The cosmic gamma-ray background rep ortedl y 

(|2008l); 



(see, for insta nce , 
Ahn fc Komatsul (l2005al lbh; 



Inoue. Totani. fc Y. Uedal 



Rasera et al.l 



2006); 



Strigari et al. (|2005l )) has a peculiarity in the en- 
ergy range 0.5 - 20 MeV. A ledge-like feature is 
visible in th e extragalactic gamma-ray s pectru m 
(Fig. 3 in IStrong. Moskalenko, fc Reimerl i|2004h ). 



The photon index here is markedly distinct from 
those of the so fte r or h arder parts of the spectrum 
llSreekumar et al] Il998l; IStrong. Moskalenko. fc Reimerl 
|2004| ; IWeidenspointner et al.l l2000f ). indicating its different 
origin. 

Moreover, this spectral band can be formed nei- 
ther by too soft emission of normal active galactic nu- 
clei, nor by too hard bla z ar-typ e AGNs contribution (see 
llnoue. Totani. fc Y. Uedal (|2008l ) and references therein). 
Attempts to consider the nuclear-decay gamma rays from 
Type la supernovae as the source have not been success- 
ful: the flux expected from the supernovae is several times 
weaker than the observed llAhn. Komatsu. fc Hoflich: [20051 : 
iRasera et al] 120061 : IStrigari et al] 120051 ). It might be well 
to point out that the precise determination of the excess 
boundaries and intensity is mod el-dependent, and the lit- 
erature values vary considerab ly (llnoue. Totani. fc Y. Uedal 
120081 : lAhn fc Komatsul l2005al fbl: IRasera et al] 120061 ). In any 
case, however, the excess becomes apparent near 0.5 MeV 
and disappears at the energies > 20 MeV. One can see that 
the energy range of the feature corresponds closely to the in- 
terval characteristic for the relic gamma emission from the 
WIMPs annihilation. This coincidence looks promising when 
it is considered that the WIMP is now one of the most prob- 
able dark matter candidates. 

At the same time, the relic gamma emission predicted 
by equations (|29[) . (|30|l . (|19|) , and (|28[) . is approximately 
five orders fainter than the observed feature (as we can 
see in Figure [2] the relic emission near the maximum at 
10 MeV has e 2 * flux ~ 10" 8 MeV/(srcm 2 s) while the to- 
tal extragalactic gamma-ray backg round at 10 MeV has 
e 2 * flux ~ 2 - 10~ 3 MeV/(srcm 2 s) (jSreekumar et al.lll998l : 
IRasera et afll2006h ). This discrepancy might result from in- 
applicability of the assumption of homogeneous dark matter 
distribution. In fact it cannot be so. The modern structure 
of the Universe appeared from some initial perturbations 
that had already existed, beyond any doubt, in the epoch 
z ~ 3 00. According to WMAP measurements (|Spergel et al] 
2003), in the recombination epoch z ~ 1100 -f- 1400 relative 
variations of the baryonic matter density were of the order 
of 10 -5 . Dark matter perturbations could be much more in- 
tensive, they were not suppressed by the radiation pressure 
in the pre-recombination epoch. Moreover, they must have 
been significantly stronger (not l ess than 10~ 3 ) to explai n 
the modern Universe structure (|Gurevich fc Zvbinl Il99rj ). 
Since the recombination happened in the matter-dominated 
epoch, the perturbations rapidly grew and at the moment 
z ~ 300 could attain very big amplitude. The presence of 
density inhomogeneities does not affect the spectrum of the 
annihilation signal but increases its intensity. This effect is 
usually described by the quantity 



C = 



<P*> 
< p > 2 



that appears as a mul tiplier in the expression f or the inten- 
sity (see, for instance. lAhn fc Komatsul (|2005al )). Of course, 
C is a function of z. 

This brings up two points: first, is it possible that 
the structure growth in the early universe proceeded so 
fast that C was as large as 10 5 — 10 6 by the moment 
z — 300? Second, if an intensive structure formation took 
place at some moment z < 300, the coefficient C(z) could 
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grow so rapidly that it far outweighed the signal diminu- 
tion owing to average density decrease. As a result a hard 
tail or even a secondary hard maximum on the energy (2- 
4)/z GeV can appear in the spectrum of the relic emission, 
which is not observed. Unfortunately, the theory of evo- 
lution of dark matter perturbations is still far from accu- 
rate. They evolved from some primordial fluctuations ex- 
isted in the very early universe. While the universe was 
radiation-dominated their growth was slow. The smallest 
perturbations were destroyed by free-streaming (for the in- 
stance of neutralinos of the m ass ~ 100 GeV this limit is es- 
timated as 10~ 12 - 10~ 6 M^ llZvbin Vvsotskv. fe Gurevichl 
ll999l ; lHofmann. Schwarz. fc Stockerll200ll )). When the uni- 
verse transits into the matter-dominated stage, the per- 
turbations start to grow rapidly, and eventually they be- 
come nonlinear and collapse. As this takes place, the small- 
est clumps collapse the first (at a time we denote as z) 
l|Gurevich fc Zvbinl [l995h . The overwhelming majority of 
these small clumps originated at that moment were subse- 
quently destroyed by tidal interaction with larger clumps 
originated later. But up to now it is the small clumps 
formed at the epoch z that makes the main contribu- 
tion to the ann ihilation rate of the dark matter since they 
are the densest dBerezinskv. Dokuchaev.~fc~ Eroshenkdl2003l ; 



iGurevich. Zvbin fc Sirotall 19971 ). At the moment z function 
C(z) underwent a rapid increase from a value close to 1 to a 
very big value. In what follows creation of the larger struc- 
tures was accompanied by smaller clump destruction, and 
C'(z) changes more smoothly. In order to be allowed to sug- 
gest that the extragalactic gamma-ray 0.5 - 20 MeV excess is 
related to the neutralino annihilation we must assume that 
z > 300, i.e. the first structures started to form before the 
universe became transparent for the photons produced by 
the annihilation. 

How realistic is such an assumption? Unfortu- 
nately, present estimations are very vague. Even the 
minimal possible clump mass for the neutralino dark 
matter is determined extremely uncert ainly (from 
10~ 1 2 M^ (IZvbin. Vvsotskv. fc Gurevichl 1 19991) to 1O _6 M 



of primordial fluctuations is not exactly flat, and the in- 
tensity of the fluctuations slightly builds up as their scale 
decreases. In the matter-dominated stage the perturbations 
w as Sp/p oc t 2 ^ 3 and a(t) o c t 2 ^ 3 , therefore Sp/p oc a 
Gurevich. Zvbin fc Sirotal 1 19971 ). If we assume that the 
small-scale perturbations left the radiation-dominated stage 
with amplitudes more than expected from the Harrison- 
Zeldovich spectrum, they collapsed earlier (in so doing we 
imply that the amplitudes of large-scale perturbations are 
fixed in such a way as they reproduce the observed large- 
scale structure of the Universe). Considering that the largest 
and the smallest clump mass scales differ by more than 20 
orders, even a small tilt of the spectrum of primordial fluctu- 
ations with respect to the Harrison-Zeldovich shape can be 
su fficient. A kindred scenario was c onsidered, for instance, 
by IGurevich. Zvbin fc Sirotal l|l997l ). In the context of the 
model examined by the authors small clumps are extremely 
dense and collapse just after the universe transition to the 
matter-dominated epoch (z > 1500). It is worthy of note 
that the situation when z > 300 can appear in much softer 
scenar ios than those similar to IGurevich. Zvbin fc Sirotal 
|l997l ). 

Let us give a more specific form to the above reasoning. 
The spectrum of primordial fluctuations is usually consid- 
ered to have a power-law shape: 



|,5 2 (k)|~fc" 



(31) 



The case when n = 1 corresponds to the Harrison-Zeldovich 
spectrum. The moment when a perturbation mode becomes 
nonli near and collapses is determined by the spectrum func- 
tion (|Gurevich. Zvbin fc Sirotalll997l ) 



T(fe) oc 



fc 3 |5 2 (k)| 



(32) 



i + {k/k eq y 

(in IGurevich. Zvbin fc Sirotal (|l997l ) F(k) is symbolized by 
F(k)). For all the perturbations we consider k k eq , and 

l^ 2 (k)| 



T(k) oc 



oc fc" 



|Hofmann. Schwarz. fc Stocked l200ll ). to say nothing of The mom ent a* = ^ when clumps collapse is 



the clumps density profile and the moment when the 
fluctuations become nonlinear. Experimental data as well 
as numerical simulations essentially cover the range of very 
large structures 10 15 - 1O 6 M (for instance, WMAP can 
observe only the biggest perturbations with masses corre- 
spon ding to a cluster of galaxies ~ 1O 15 M0 (|Spergel et al.l 
120031 )). Properties of smaller clumps are us ually obtained 
by approximation (|Ahn fc Komatsul l2005al ). However, in 
order to obtain any parameters for the tiny clumps of mass 
10" 6 - 1O" 12 M one has to extend the results by 12-18 
orders. Another source of uncertainties is the spectrum of 
primordial fluctuations. Usually it is deemed that it has 
flat Harrison-Zeldovich shape. In this case the moment of 
the first intensive clump creation is estimated as z ~ 80, 
though some individual clumps c ollapsed much earlier 
|Green. Hofmann. fc Schwarzl 120051 ). If such a scenario 
was indeed realized, the interpretation of 0.5 - 20 MeV 
excess as a result of neutralino annihilation is out of 
the question: otherwise a strong maximum at the energy 
4/.S GeV ~ 50 MeV would appear (which contradicts to the 
observations) . 

On the other hand, let us suppose that the spectrum 



oc k 



(33) 



(34) 



which correlates with the results of iBullock et al.1 (|200ll ). 
In the case of Harrison-Zeldovich spectrum F(k) is flat; the 
structures of various scales appear almost simultaneously 
and have similar densities. If n > 1, smaller clumps appear 
much earlier and have much higher density (since the uni- 
verse density is much higher at the moment they collapse). 

We can now construct a toy model of structure forma- 
tion with the help of an ap proach used in (|Ahn fc Komatsul 
l2005al : lBullock et al.ll200lT ). We take |c5 2 (k)| ~ k 2 , i.e. n = 2, 
and the minimal possible clump mass M m i„ — 1O _7 M0. 
Since the clamp mass relating to a perturbation mode is 
proportional to k~ 3 , we obtain from (|34p : 



M 



M 
M 



(35) 



Here we have introduced the maximum perturbation mass 
M ma x that collaps e s at a = 1. This equation is consistent 
with IBullock et all i|200ll ). where it was accepted M max ~ 
1.5 • 10 13 /i _1 M Q (ft = 0.7). Then the first clumps in our 
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model appear at 2 = i.e z — 2450. We can rewrite (f34 



M 

7 



(36) 



In accordance with lAhn fc Komatsul (|2005al ) we can repre- 
sent the boost factor C(z) as a product of three multipliers: 



C(z) = A(z)-F coll (z)-[C hal °] 



(37) 



(see all the details of the model and the notation in 
(|Ahn fc Komatsul l2005al ; iBullock et ail 120011 )1. For a flat 
ACDM universe A(z) = (18tt 2 + 82a; - 39x 2 )/(x + 1), where 
x — (pm/ per) — 1 and p m , p cr are the matter and the critical 
universe densities at given z. 

The matter fraction collapsed into cosmological halos 

is: 



(38) 



Hereafter M* is the maximum having col- 

lapsed at given z^ The oretic al model 

ijBerezinskv. Dokuchaev. fc Eroshenkol 120031') an d nu - 
merical simulations (jDiemand. Moore fc Stadell 120051 ) 
give for the differential number density of clumps in the 
comoving frame of reference -j^ oc M~ 2 . Substituting it in 
(|38[) and using ((36]), we obtain: 



F u 



F coU (z) = -^ln , 

46.8 VM m „, 



M* 



7.8 



(39) 



wh ere F^ nll = F co g (a = 11 We will adopt the value obtained 
bv lAhn fc Komatsul (|2005al ) F° oll ~ 0.8. 

The factor [C hal °] = / C halo Mj^dM/ j dMj^M rep- 
resents the "halo clumping". If we adopt the Navarro- Frenk- 
White clump profile than 



cg, r (l-l/(l + c mr ) 3 ) 

9(ln(l + C vir ) — Cvir /{I +Cvir)) 2 



(40) 



In order to describe the halo concentratio n parameter c v j r 
evoluti o n we will use an eq uation from lAhn fc Komatsul 
j2005al ): lBullock et all (|200ll ) with K = 8, F = 0.01 



c vir {M,a) = Ka ( F— 



= K~ , 

a V FM 



(41) 



This formula was obtained as a fit of N-body simulations 
and is valid only for a limited range of M and z covered by 
them. We have to consider a much wider mass and red shift 
range, and (|4ip gives too large c v i r for the smallest clumps 
(for instance, for the minimal mass clumps at the present 
epoch c V i r — 20000). Such a huge value seems unlikely and 
indicates that eq u ation (1411) should be corrected. Following 
lAhn fc Komatsul (|2005al ). we assume that for any chosen 
clump mass the concentration parameter rises only up to 
Cvir = 100, and after remains constant. 

Now the toy model is defined. The boost factor C(z) 
evolution curve predicted by it is represented in Figure [3] 
We can make two principal conclusions. First of all, z = 2450 
and C(z — 300) ~ 1.5 ■ 10 5 in the model considered, i.e. the 
first structures appear very early (z 3> 300), and at the mo- 
ment z = 300 the boost factor is large enough to explain the 
discrepancy between the observed and the predicted signal 
intensities. Second, the boost factor C(z) grows respectively 



o 




Figure 3. The boost factor evolution with z. 



slowly after z = 300. From the moment when the universe 
became transparent for the annihilation photons up to the 
present moment it increases only on an order, while the sig- 
nal without regard for the dumpiness rapidly falls as z 2 ' 5 
owing to the dark matter density decrease ((9}. It means 
that almost all the signal appears at the moment z ~ 300 
and the resulting spectrum has neither a hard tail nor a 
secondary harder maximum. In closing we remark that the 
foregoing structure formation scheme is no more than a toy 
model presented here only to illustrate that the structures 
in the universe could appear early, and the interpretation of 
the excess as a result of WIMP annihilation is principally 
possible. 

The dark matter annihilation has already been in- 
voked to explain the 0.5 - 2 MeV excess (|Ahn fc Komatsul 
l2005al lbl: iRasera et al.ll2006f ). Since the photon spectrum is 
relatively soft, the authors introduced a low-mass dark mat- 
ter candidate (Mdmp < 100 MeV). Its annihilation cross- 
section is sizable t o provide the observed signal (crv) ~ 
2.5 • 10~ 26 (cm 3 /s) (|Ahn fc KomatsvJl2005al ) that is at least 
no less than the typical weak interaction cross-section at this 
energy scale. Such a low-mass dark matter candidate with 
such a significant cross-section is now ruled out (by the ac- 
celerator experiments) in ordinary schemes like MSSM, but 
there is an interesting possibility to introduce it in more 
sophi sticated scenarios (jBoehm et al 1 l2004l : iBoehm fc Favetl 
|2004| ). 

In our case, we can manage with usual heavy WIMP 
candidates. As we could see, the relic gamma-ray signal red- 
shift is z ~ 300. Originally, the photons had the energy 
1 -T- 5 GeV and they were produced by the annihilation of 
ordinary WIMP particles like the lightest neutralinos. On 
the other hand, if our interpretation is true, it counts in 
favour of a relatively light WIMP (Mdmp ~ 100 GeV). If 
Mdmp S> 100 GeV, the typical energy of the producing 
photons is higher, and the feature in the spectrum must be 
harder than the observed one. 

Thus, the critical point for the 0.5 - 20 MeV excess 
interpretation as a result of relic neutralino annihilation is 
determination of the moment z when the first structures 
in the universe appeared. If the structures had appeared 
before the universe became transparent for the annihilation 
products (i.e. z > 300), then there are strong arguments to 
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believe that the excess 0.5 - 20 MeV is created by the relic 
WIMPs annihilation. Above all, the characteristic energies 
of the spectra agree. Besides, the WIMP is now one of the 
most probable dark matter candidates, and the coincidence 
does not look random. The discrepancy of the predicted and 
observed signal intensities can be naturally explained by the 
nonuniform structure of the dark matter. If the first clumps 
appeared later (z < 300), the excess undoubtedly could not 
be produced by neutralino annihilation. Further progress in 
the Universe structure formation understanding will be able 
to shed light on this problem. 
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